Spatial Entanglement of a Free Bosonic Field 
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In this paper we discuss the entanglement properties of a thermal non-relativistic free bosonic 
field. We demonstrate how to formally construct spatial modes in order to use a continuous variable 
separability criterion and show that the construction of the modes plays a significant role in the way 
the entanglement manifests itself. For instance the presence of entanglement depends on the shape 
of the modes and their separation. The temperature of the gas is another crucial factor and for one 
choice of modes we show that entanglement can be found at arbitrarily high temperatures as long 
as we can zoom in on suitably small regions in space. Moreover, we show that the entanglement 
here is useful as it can be extracted to a pair of localised systems. 



I. INTRODUCTION 

Entanglement, a fundamental ingredient of quantum 
mechanics, leads to correlations between subsystems that 
are greater than anything possible classically and is the 
main resource for quantum computations. However these 
quantum correlations are delicate and can easily be de- 
stroyed. Therefore it is logical to hunt out systems that 
possess an amount of natural thermal entanglement when 
they are under certain conditions, for instance below a 
certain temperature, so that if one keeps the conditions 
stable, the entanglement definitely exists even though the 
system may be open. 

Unlike previous papers that use macroscopic variables 
to witness entanglement in discrete systems - for instance 
temperature [|| , heat capacity Q and magnetic suscepti- 
bility we will use an existing separability criterion 
to investigate whether a Bose field is entangled and look 
at how such variables affect the entanglement. This is 
useful to know as entanglement found at higher temper- 
atures is easier to obtain than in situations where the 
exact ground state is vital 0, [l(| . 

We will focus on spatial entanglement, which means 
that localised regions of a system exhibit particle num- 
ber correlations that arc non-local. When considering 
spatial entanglement it is necessary to work in the second 
quantised formalism, as entanglement in many-body sys- 
tems can exist between the occupation number of modes 
and not just between internal degrees of freedom of the 
particles themselves. This allows for the interesting case 
of single particle entanglement, which was elsewhere (2C| 
shown to be of the same form as the usual EPR entangle- 
ment. Crucially, the first quantisation is misleading since 
it can generate inconsistent conclusions. If we worked in 
the first quantised picture, defined two spatially sepa- 
rated regions, the Hilbert space would not be a direct 
product of the Hilbert spaces of the two regions and can- 
not be decomposed into spatial modes. Thus we must 
work in the second quantised picture where each region 
represents an independent mode and our Hilbert space 
is a tensor product of the two. We will proceed to in- 
vestigate entanglement of the free bosonic field in this 



manner. This is an interesting problem as not only will 
we show that spatial entanglement can be extracted for 
computational purposes but the physics here is curious as 
we find an underlying entangled structure in a relatively 
simple system. 

To motivate our work further note that a recent paper 
pd| provided a general argument for why energy might 
act as a good witness for spatial entanglement in CV 
many-body systems and demonstrated that the critical 
temperature for this entanglement is remarkably similar 
to the critical temperature for Bosc-Einstcin condensa- 
tion. In the present paper we put foundations to this 
problem by tackling it from an entirely different direc- 
tion. Wc formally construct spatial modes in terms of 
ladder operators for regions of space and apply an exist- 
ing separability criterion [l3| which unambiguously de- 
cides whether the two spatial modes are entangled. Un- 
like [Til ], in this constructive scenario, it is clear how the 
entanglement emerges and our criterion is both necessary 
and sufficient for bi-partite divisions in its most general 
form. Additionally, wc will discuss the results in terms of 
thermal states, give an extraction procedure and contrast 
our methods to another related paper [ljj ■ 

The paper is organised as follows. Let us proceed by 
firstly providing an overview of Gaussian states and the 
separability criterion that we will use to detect the en- 
tanglement in section [ill Therefore if the reader is fa- 
miliar with Gaussian states please go directly to section 
Mil We next analyse the spatial entanglement in the gas 
and discuss the role of the detector profile in section IIIII 
Finally we show that the entanglement is useful as it can 
be extracted to a pair of localised systems in section IIVI 



II. GAUSSIAN STATES 

Gaussian states are common in the physical arena as 
they closely approximate coherent and squeezed states of 
light fl7| and thermal states of a quadratic Hamiltonian 
are Gaussian. The set of Gaussian states have Gaussian 
characteristic functions and are therefore expressed en- 
tirely in terms of their first and second moments. We 
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will now develop a mathematical formalism to allow us 
to look at two-mode Gaussian states closer. Please note 
that complete reviews of Gaussian states can be found in 

A quantum bosonic field is equivalent to a (possibly in- 
finite) set of independent quantum harmonic oscillators. 
Thus for a known number of bosons the Hilbert space 
H = (£) 2 k=1 H-k of such a system is a tensor product of 
single-oscillator Fock spaces Ttk- The creation and anni- 
hilation operators 4 and dk acting on the Fock spaces 
can then be related to the phase-space operators through 



2mujk 



(4 + a k ), 



Pk 



hmLUk 



(4 " & k), 



(1) 

where ujk is the angular frequency of the fcth mode. If 
the operators are arranged in a vector as follows 



(2) 



the canonical commutation relations can be expressed in 
the compact form as 



(3) 



where the subscripts i,j = 1,2,3,4 refer to the different 
components of this vector. The matrix f2y is defined as 



J = 



k=l 



1 

-1 



(4) 



It is well known that the state of the gas is described 
by the positive trace-class density operator p that is used 
to calculate expectation values of the phase-space opera- 
tors, which form the first and second moments describing 
the Gaussian state. These can be grouped together in the 
covariance matrix (CM), which is defined through its el- 
ements as 



(5) 



where (O) = Tr[Op]. The CM 7 is real, symmetric and 
holds all the available information about the state of the 
gas and like the density operator p the CM must fulfill 
certain criteria to be a bona- fide CM. Note that the dis- 
placements can be altered by local operations and 
classical communications (LOCC) and do not affect the 
entanglement properties of the system and are set to zero 
here without loss of generality. 

Thus, using © and §5§ is it straight forward to show 
that the CM must obey following inequality 



7' 



in > 0. 



(0) 



This inequality is the only requirement for the CM of the 
gas to represent a physical state. 

It will be useful to manipulate the CM in order to 
simplify the calculations. Using local symplectic trans- 
formations Si and S*2 such that S = Si © S2 preserve 



any entanglement properties of the CM, a ten parameter 
CM of two modes can be reduced to only four param- 
eters. The CM is now in its standard form 7 SF and is 
expressed as 



S T -/S = 7 



SF 



fa c 

a d 

c b 

\0 d b 



(7) 



There are four symplectic invariants associated with the 
CM, namely det7 = {ab-c 2 ){ab-d 2 ), detA = a 2 ,detB = 
b 2 and detC = cd. These invariants are unique up to the 
relative sign of c and d for any covariance matrix 7. The 
symplectic eigenvalues EH EH arc defined as the ones 
for 7 itself and are expressed in terms of the invariant 
quantities (a, b, ...etc.) as 



v± 



1 
71 



+ b 2 + 2cd ± y/(a 2 + b 2 + 2cd) 2 - 4det 7 . 

(8) 

The uncertainty relations on the CM are satisfied if v_ > 
1. It is known that only pure states can saturate this 
bound while mixed states necessarily remain above it. 

When dealing with entanglement it is often useful to 
calculate the purity of the state under consideration be- 
cause as a states purity increases its entanglement usually 
increases too . When using a density operator p we can 
calculate the purity through p = Tr[p 2 ], which translates 
into 



M = Tr[p 2 ] = 



1 



Vdet 7 



(9) 



where for pure states det 7 = 1 and det 7 > 1 for mixed 
states. 

To investigate entanglement of a two-mode Gaussian 
state we use a criterion by Simon [l3| that defines the par- 
tial transpose operation for CV states. The transposition 
operator used to detect entanglement in El, 0| is equiva- 
lent to partial complex conjugation (or time reversal) of 
a CV system. Simon saw this as a 'mirror reflection in 
phase space', however, a mirror reflection in phase space 
is not a canonical transformation and cannot be physi- 
cally implemented. It is used here as a mathematical way 
of detecting entanglement. 

To apply this operation to a bi-partite CV state we 
change the sign of the momenta of one constituent sub- 
system. As in El, [2| when the transposition is applied 
to one part of an entangled system it causes the state 
to become unphysical, the uncertainty relations ^ with 
our new CM, 7 Ta will become unphysical if the state 
is entangled. In fact, it is straightforward to check this 
relation using the symplectic eigenvalues ([5]) and i>_ > 1. 

In addition it is useful to know the amount of entangle- 
ment in a system, which can be calculated through the 
logarithmic negativity Q 



Em = max{0, - log 2 (^_)}. 



(10) 



The logarithmic negativity Ejj is therefore zero for sep- 
arable states and increases unboundedly as i/_ — > 0. 



3 




02 (x) 



<f>i(x) 

(f) (x) 



R 



AS 



Q 



L 



FIG. 1: Visualisation of the system, in which a thermal non- 
interacting bosonic gas is trapped in a ID infinite square well 
of length L and the particles occupy decoupled momentum 
modes <j>k(x). We wish to investigate entanglement between 
two spatial modes R and Q (grey shading) with a finite dis- 
tance between them, denoted by AS. 



III. ENTANGLEMENT BETWEEN SPATIAL 
MODES 

The formalism in the previous section allows one to cal- 
culate the entanglement between two regions of a field, 
but first, the regions need to be constructed so that the 
separability criterion from the previous section is appli- 
cable. 



A. Defining Spatial Modes 

Before we construct the spatial modes (regions) them- 
selves we shall briefly describe in more detail the sys- 
tem under consideration. For simplicity we take a ID 
infinite square well potential of length L, see Fig. [TJ 
filled with a free bosonic gas in a thermal state p = 
cxp(— (H — fj,N) /ksT) I 'Z at temperature T. The quan- 
tities H, ks, fJ- and N are the Hamiltonian, Boltzmann's 
constant, the chemical potential and the number of par- 
ticles in the system, respectively. The grand partition 
function is Z = ir[exp(— j3(H — /J.N)}. Note, that confin- 
ing geometries of higher dimensions could be considered 
but this does not contirbute much to the present evalua- 
tion. In momentum representation, the Hamiltonian of a 
non-interacting bosonic gas in its second quantised form 
is simply H = VJ fc Eka\a,k, where k are the momentum 
modes of the system , and a k are the annihilation 
and creation operators of particles in these modes and 
Ek = fiujk with the usual dispersion relation for massive 
particles u>k = hk 2 /2m. As the Hamiltonian is quadratic 
the thermal state is Gaussian and can be described en- 
tirely by its first and second moments. 

The Hamiltonian describes the occupation of decou- 
pled momentum modes (j>k(%) = \f^jL sin(fcx) with 
k = nl/ L and consequently the state p is separable w.r.t. 



a decomposition into momentum modes. Nevertheless 
entanglement could exist in position space and we wish 
to investigate any spatial correlations. To apply the CV 
separability criterion (l3| . position and momentum op- 
erators for the regions should be constructed. In lattice 
systems the idea of spatial modes is clearly defined in 
terms of individual or groups of lattice sites and the cor- 
responding operators would be the position and momen- 
tum of the particles at the lattice sites. However, in a 
CV system one has the freedom to include any region of 
space in the spatial mode, which may in turn contain any 
number of particles. The modes are then constructed by 
averaging position and momentum operators with a com- 
plex detector profile over the set of points belonging to 
the region. Although our operators must be position and 
momentum, they will therefore be the averaged position 
and momentum of all the particles (or excitations) in the 
spatial modes and not the position and momentum of the 
particles themselves. 

In fact, position and momentum operators acting on 
Fock spaces of individual momentum modes were defined 
in the previous section ([T]). So, in order to construct the 
spatial modes here we firstly take those phase-space oper- 
ators ([T]) and transform into position representation with 
the discrete Fourier transform u(x) = Y]u—\ (j^k{x) Uk and 
similar for p{x) . And secondly, as space - unlike momen- 
tum in our system - is continuous and a point in space 
is of zero volume and hence meaningless to talk about, 
we investigate entanglement between spatial regions of a 
finite size, denoted by R and Q (Fig. [1]). The conjugate 
operators u(x),p(x) become ur,pr,uq and pq by aver- 
aging over a detector profile g(x) localised in R and Q, 
for example the operator ur is ur = J R dx g{x)u{x). 

We can take any normalised complex function for g{x) 
and of course there are an infinite number to choose from. 
If we find just one g{x) for which the separability con- 
dition is violated, the gas is spatially entangled. On the 
other hand, the gas is only truly separable if there is 
no detector profile that gives entanglement, but to con- 
firm this one would have to search over infinitely many 
g(x). Naturally our operators must obey the usual com- 
mutation relations [ur,Pq] = [ur,uq] = [Pr,Pq] = 
and [ur,pr] = ih, which determine the normalisation of 
9(x). 

Now that we have constructed the operators that de- 
fine our spatial regions, we can apply the separability 
criterion [l3| to the gas. 



B. Applying The Separability Criterion 

It is straight-forward to check that the transformation 
between the momentum and spatial modes is linear and 
the Gaussian nature of the state in the new canonical 
operators is preserved. This means that the state can 
be described fully by a two-mode CM ([5]), with the two 
modes representing regions R and Q. 

By expanding the operators in terms of the recipro- 



4 



cal space annihilation and creation operators, and a\, 
the constituents of the covariance matrix can be evalu- 
ated easily in the momentum Fock basis. The resulting 
CM is of the form 



/ A E 

B F 

E C 

V F D 



(11) 



For completeness the general form of the CM elements 



A,C= 2(u 2 R ) f 



■[(/, 



2tr|( / dxg(x)u(x)) (12) 
fc 



B,D = 2(f RiQ ) p = ^7r 2 /c 2 (Z^ Q ) 2 coth 



e = ({u r ,u q }) = ^;_L_i*i* 



coth 



(£fc - M)/3 



F={{p R ,p Q }) - ^7r 2 fc 2 X^coth 



where X|j 



= J R q dx g(x)ipk(x). The factor ^/h/2mu>k 
in the operator results in the l/ir 2 k 2 in A, C and 
_E above ((T^J) and terms containing fc 2 occur in the 
sums for the momenta uncertainties (elements B, D and 
F). These terms play an important role in the conver- 
gence/divergence properties of the elements and conse- 
quently in the resulting entanglement. 

Everything is now expressed suitably for the CV sep- 
arability criterion described in section ITT] We would like 
to arrive at an inequality that must be positive for sep- 
arable states and negative for entangled states so let us 
apply the partial transposition which means F — > —F. 
Examining the symplectic eigenvalue of the partially 
transposed CM leads to the following expression in terms 
of the original CM elements, 

£ = 1 + (AC - F 2 ){BD - E 2 ) -AB-CD + 2EF. (13) 

If £ > the two spatial modes are separable, on the other 
hand £ < means that the modes are entangled. 

However, before we can even speak about entangle- 
ment between the spatial modes the CM must represent 
a physical state, which is the case when the uncertainty 
relation (J6]) for the original CM 7 is satisfied. 



C. Choosing A Detector Profile 

Before we can check the entanglement of the gas 
through (fT3|) , a specific form for the so-called detector 
profile g(x) must be assumed. In other words how should 
the set of points belonging to the regions R and Q be av- 
eraged over? 



1. Top Hat Detector Profile 

Let us start simply by giving equal weight to all the 
points belonging to the regions i.e. we use a top-hat 
function that gives support inside the regions and is zero 
outside of them. This is a natural choice as in practice 
many detectors click if something is in a certain range but 
will not weight any two positions in that range differently. 

With this choice of detector profile the system is al- 
ways separable w.r.t. the two spatial modes, indepen- 
dently of the size of the regions, their separation and the 
temperature. This is intuitively clear as in a thermal 
state all possible momenta are mixed and correlations 
are lost. Mathematically this is clear as the momentum 
entries (P R ) etc in our CM contain terms proportional to 
k 2 in the sums over k and will not converge. To combat 
this and arrive at a less mixed state, the state could be 
prepared so that it contains only a particular portion of 
the mixture. However entanglement can be found with- 
out manipulating the state at all, so let us look firstly at 
this. Although later we will come back to using a top 
hat function as the detector profile and allow for some 
selection of momenta (see section UlI C 3p . because inter- 
estingly in that case entanglement is found to exist, under 
some circumstances, at very high temperatures. 



2. Gaussian Detector Profile 

In the previous section the detector profile g(x) was 
chosen to be a top hat function, however this meant that 
all momenta were detected with equal weighting and no 
entanglement was uncovered. Conversely, one could have 
chosen g(x) so that it naturally gave a lower weighting to 
higher momentum modes so that the state becomes less 
mixed. This is what we will try next. 

The freedom to choose the detector profile illustrates 
that there are infinitely many choices of modes between 
which there could be entanglement. This is because a 
free gas has infinitely many degrees of freedom. If we 
find one pair of spatial modes that are entangled than 
we can safely say that the gas is entangled. One could 
choose any detector profile, but by fixing the detector 
profile the entanglement between any other combinations 
of modes is excluded and is therefore less than the actual 
amount of entanglement in the gas. 

Naturally it is difficult to find the g(x) that picks 
up the entanglement, so instead we start by tak- 
ing the function T r q = J R q dx g(x)(j>k{x), which oc- 
curs in (fT2"|) and see that for the top hat detector 
profile it was T RQ = C/ky/2/L{cos{Ri{Qi)k/ L) - 

cos(i?2(Q2)fc/-k)); where i?i(Qi) and i?2(Q2) denote the 
edges of the regions R or Q. We use this as the 
starting point to define a new I R q and combine it 
with a function that gives greater weighting to a set of 
low momenta depending on the size of the region un- 
der consideration. Thus 1 RQ = C(cos(Ri(Qi)k/L) — 
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FIG. 2: The plot shows the logarithmic negativity for the 
region size Ax = L/W against the separation of the regions 
and the temperature of the gas. The amount of entanglement 
increases as the temperature is lowered and only exists when 
the regions are sufficiently close to one another. 

cos(R2(Q2)k / L)) exp(— k 2 Ax 2 ), where the normalisation 
C is determined through the commutation relations. 
Here, although the detector profile has not been specified 
directly, we have chosen one that influences the range of 
momenta it picks up Afc depending on the size of the 
region Ax. 

The new choice of detector profile has the desired effect 
as entanglement is now found in the gas without any 
other manipulations. 

Let us look now at how the entanglement manifests 
itself. Firstly, unlike for a top hat detector profile if mo- 
menta in the state were pre-selected (as discussed in the 
following section Ull C 3[) . here the entanglement cannot 
exist at arbitrarily high temperatures even if the regions 
are very small as in this case the uncertainty relations 
© arc not satisfied. Entanglement exists for regions in 
the size range of Ax ~ L/20 — > L/2 and the amount 
of entanglement for regions of the size Ax = L/10 and 
for non-zero temperatures of the gas is plotted in Fig. 
[5] using the logarithmic negativity as a function of their 
separation. There is no entanglement above a certain 
temperature and when the regions are far away from one 
another. 



3. Top Hat Detector Profile and Selecting Momenta 

It is clear that a thermal Bosonic field is entangled 
when a suitable detector profile that limits the mixed- 
ness of the gas is chosen. In order to compare the present 
work with and illustrate the interesting case of en- 
tanglement at very high temperatures, we would like to 
show how one can still find entanglement in the gas with 
top hat detector profile by preparing the state with the 



necessary mixedncss. 

It is straight-forward to examine the mixedness of the 
gas using the purity relation given in eq. © as the en- 
tanglement criterion (|13[) can be written in terms of the 
purity 



y/AB + CD- 2EF - 1 ' 

The gas can be prepared with a greater purity by pro- 
jecting out all configurations over a certain momenta. 
Once the purity satisfies (fT4|) the gas will be entangled. 
However, with a top hat detector profile this action can- 
not be performed locally as we cannot distinguish the 
momentum modes in the spatial regions alone as they 
are only orthonormal w.r.t. the entire well. The pro- 
jection is therefore a global operation, but this is just 
part of the state preparation as we are defining the shape 
of the spatial modes upon which the entanglement very 
much depends. Indeed, the projection proves to be a 
profitable exercise as afterwards £ < and we uncover 
spatial entanglement, under certain circumstances, be- 
tween the two spatial modes. 

We know that for entanglement to emerge, the gas has 
to be of a certain purity which means that the state must 
only include particles below a given momenta. But how 
does one estimate the maximum momenta that an entan- 
gled field should have? It turns out that the maximum 
momenta for an entangled set of particles is determined 
by Hcisenbcrg's uncertainty relations. It was noted in the 
previous section that the CM must fulfill the uncertainty 
relation ^ if it is to represent a physical state. In simple 
terms, this means that the state should occupy an area 
in phase-space greater than Ax Afc ~ 1, where Ax is the 
size of the regions and Afc is the number of different mo- 
menta in the state. When defining the spatial regions we 
have limited the size of our Ax so consequently we must 
ensure that the state is prepared with an adequate Afc 
to satisfy this relation and that we do not try to measure 
both quadratures with a greater precision than possible. 
However, only by considering states that have the small- 
est possible uncertainty is the entanglement uncovered 
(Fig. [3]). Fig. [3] indicates how many momenta have to be 
included in our state for the two regions to be entangled 
at a fixed T and confirms that Afc ~ 1/Ax is roughly 
obeyed. 

The entanglement depends upon the width of the re- 
gions, their separation and the temperature of the sys- 
tem. Entanglement for given sized regions exists when 
the regions are suitably close and below a certain critical 
temperature Tc- If we investigate the nature of this crit- 
ical temperature with respect to the size of the regions, 
from Fig. 0]it is evident that the smaller the regions, the 
higher the critical temperature for entanglement. In fact 
we assert that at very small length scales entanglement 
exists at arbitrarily high temperatures as the critical tem- 
perature behaves roughly as Tc oc Ai"i with the size 
of the regions Ax. The high temperature entanglement 
occurs between two very small regions with a top hat 
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FIG. 3: We plot the number of momenta (Ak) that is included 
in an entangled state against the size of the regions Axr = 
Axq = Ax when the two regions are touching at L/2. We 
sum from the lowest momentum k — n/L up to k m i n where 
one first finds entangled states (this is our AA; mm ). For some 
Ax, entangled states include a wider range of momenta from 
k = n/L up to kmax (Afcmax). It is evident that for smaller 
regions entangled states must include an increasing number 
of modes in order to satisfy the uncertainty relations. The 
shape of this graph illustrates that Ak ~ 1/ Ax is obeyed. 



detector profile precisely because of the uncertainty re- 
lations (J5]). On one hand, the Gaussian detector profile 
suppressed the occupation of higher momentum modes at 
high temperatures which meant the uncertainty relations 
were not satisfied. Whereas here it is always possible to 
satisfy the uncertainty relations, but in order to generate 
a Afc large enough to do so, the gas has to be heated up 
to a finite T, which results in high temperature entangle- 
ment. 

Regardless of the choice of detector profile any spa- 
tial entanglement in the gas is in the form of particle 
number correlations between the two spatial modes and 
as it is not between the particles themselves, one may 
question whether the entanglement is genuine and use- 
ful or whether it is created by artificially defining the 
regions. The following section demonstrates that this is 
genuine entanglement by describing an extraction process 
whereby the entanglement is transfered to two localised 
systems which could then be used for teleportation. 



IV. EXTRACTING ENTANGLEMENT FROM A 
BOSONIC GAS 

We know that the entanglement in the gas is strongly 
dependent on the choice of detector profile g(x) that de- 
fines the shape of our region. However, in this section the 
detector profiles will result from the effective volume of 
two localised systems (here we take atoms but we could 
of course take quantum dots or any other viable system), 
which will interact locally with the gas. Entanglement 
will be extracted from the gas so not only will we have 
proper valuable EPR-type entanglement, but the extrac- 



FIG. 4: The critical temperature Tc as a function of the 
size of the regions, where Axr — Axq and the regions are 
touching at L/2. For larger regions the critical temperature 
is very low, but for smaller and smaller regions Tc increases. 



tion also provides another way of locally reducing the 
number of momenta in our state. 

The method we use for extraction is similar to a gen- 
eralised scheme by Kaszlikowski et al. [24| although here 
we use a different Hamiltonian. Our set up is as follows. 
Two systems - held by Alice (Sa) and Bob (Sb) - are 
placed in close proximity to the gas and their natural fi- 
nite width is used to define the shape (detector profile 
g(x)) of two separate spatial regions. The systems are 
localised in two distinct regions of space and start in a 
separable state with respect to one another. The probes 
will interact with the gas with a sufficiently short time 
in order to deny them from sharing information with one 
another via the bosonic field. 

The Hamiltonian that provides this interaction is 
H(t) = T(t)(u R Ps A + uqPs b ), where T(t) is the interac- 
tion strength. The position of the gas located in region R 
couples uniquely with the momentum of Alice's system 
Ps A and the position of the gas in the other region, Q, 
couples locally to the momentum of Bob's system Ps B . 
To select the entangled set of momenta the probe atoms 
should have a gaussian shape to them which mimics the 
g(x) used in the previous section. This would be the case 
for instance if an atom was in the lowest energy level in 
a harmonic trap. 

To show that the probe systems are entangled after we 
have applied the interaction for some time t, we must 
determine the time evolution of the state. The sys- 
tem is initially in the state p(Q) = p ff |00)(00| (where 
1 00) = |0)a <& |0}b and here p g is the density opera- 
tor of the bosonic gas) and as the full time evolution 
is difficult to establish, as we have chosen a suitably 
short interaction, we can calculate the state at time 
St, perturbatively, to first order. The state of the gas 
and probes as limat^or = / ' dtT(t) « 1 is p(5t — > 
0) = p(0)+T 2 hmu;/2(u R p g u R \10){10\+u Q p g u Q \01)(01\ + 
(u R p g u Q \10){01\+h.c.)) +T 4 u R u Q p g u R u Q \ll) (11|. 

We would like to know whether the probe systems have 
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now become entangled in momentum and for this we 
trace out the gas and compute the Peres-Horodecki cri- 
terion 0, 0]. Before we apply the partial transposition 
the density matrix of the atoms, written in the basis of 
the lowest two momentum states, is of the form 

/ 1 \ 
_ I x y I 

P Atoms - y z , 

\0 5/ 

where x oc tx\puj^,z oc tr[pug], y oc tr[puRUQ] and the 
term S oc tr [pu^Uq] is very small (as it its proportional 
to r 4 ) but non-negligible. If y > 1/2^5(5 + 1) is satis- 
fied the atoms will become entangled and as i5 < 1 wc 
can approximate this by y > 0. If the off-diagonal ele- 
ment y exists then we can extract entanglement from the 
gas. The term y octr^UflUg] is the overlap between the 
positions of regions R and Q and its existence depends 
upon g(x) the temperature of the gas and the separation 
of the two regions. 

V. DISCUSSION AND CONCLUSIONS 

Let us now compare our findings to two other papers 
that investigated field theoretic entanglement. 

Firstly we found that entangled fields had a maximum 
energy. If we compare this to the results of [ll[ , it is not 
surprising that the minimum energy of the separable set 
of states is fixed by the uncertainty relation. There, the 
authors used the uncertainty relations applied to regions 
of space to find the minimum energy that separable states 
can have and said that states with a lower energy are en- 
tangled. Unlike [ll[, here we have tackled the problem 
constructively by defining proper spatial modes and ap- 
plying a well known separability criterion and as a result 
we can see clearly that by restricting ourselves to slow 
moving particles in the two regions the state becomes 
purer and therefore becomes entangled. 

Secondly, another paper [l2j that analyses entangle- 
ment of a field, constructs spatial modes quite differently 
to us. In [l2| each region of space was mathematically 
converted to a two level system and the field was pro- 
jected onto these two states thereby reducing an infinite 



dimension system down to just two degrees of freedom, 
which allowed the application of the Peres-Horodecki cri- 
terion [l], 0] • On the other hand we reduce our field to a 
harmonic oscillator for each region and allow each one to 
be populated by bosons. The latter is more operational 
than |l2| as the way we average over the set of points in 
each region (i.e. we must specify a detector profile for 
each region) proves to be crucial for revealing entangle- 
ment. 

In this paper we have shown that the remarkably 
simple system of non-interacting bosons is entangled in 
space. Unlike discrete systems such as spin chains that 
have their spatial modes naturally defined here one must 
create spatial modes by averaging operators over a set 
of points with some detector profile g(x). It is there- 
fore natural that the choice of g(x) effects the entangle- 
ment that is uncovered. Thus the entanglement resulting 
from two different choices of spatial modes has been anal- 
ysed. Firstly, if one averages over all points in the regions 
equally (i.e. the detector profile is a top hat function), 
to find entanglement one has to prepare the gas by low- 
ering the number of momenta in the state using a global 
operation. In this case however, entanglement exists at 
very high temperatures if we consider very small regions 
due to the uncertainty relations. Secondly, a Gaussian 
type detector profile that weighted the higher momentum 
modes less than the lower ones uncovered entanglement 
locally. In this case, the entanglement always existed be- 
low a certain temperature and disappeared altogether for 
very small regions. Moreover we have demonstrated that 
the entanglement is useful as it can be extracted from 
the gas. 
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